function [dfitC,flux00,fluxAll,diffF]=predictMIRNA_FC_single_flux(modelCancer, CancerLiver_exp, miRNAtarget, ParsedGPR, Prxn)

N=size(miRNAtarget,1);

dfitC=zeros(N,1);
  

opFC=optimizeCbModel(modelCancer); flux00=opFC.x;


i=0;
t=cputime;

for i=1:N
    i
    
    for j=1:size(miRNAtarget,2)
        if isnan(miRNAtarget(i,j))
            break;
        else
            temptarget(j,1)=miRNAtarget(i,j);
        end
    end
    
    mi_exp_geneLocus=temptarget;
    mi_exp_geneFDR=ones(size(temptarget));
    mi_exp_geneFC=ones(size(temptarget))*(-1);
    
    
  
    
    [miCancer]=UB_update_withFC(modelCancer,mi_exp_geneFC,mi_exp_geneFDR,mi_exp_geneLocus,CancerLiver_exp,ParsedGPR,Prxn);

    

    
    opFCm=optimizeCbModel(miCancer);fluxAll(:,i)=opFCm.x;diffF(:,i)=opFC.x-opFCm.x;
    
    
    dfitC(i,1)=opFC.f-opFCm.f;
  
    

    
end
cputime-t